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Abstract. 

We present our recently developed 3-dlmensional chemodynamical code for galaxy 
evolution. This code follows the evolution of different galactic components like stars, 
dark matter and different components of the interstellar medium (ISM), i.e. a diffuse 
gaseous phase and the molecular clouds. Stars and dark matter are treated as colli- 
sionless N-body systems. The ISM is numerically described by a smoothed particle 
hydrodynamics (SPH) approach for the diffuse gas and a sticky particle scheme for 
the molecular clouds. Additionally, the galactic components are coupled by several 
phase transitions like star formation, stellar death or condensation and evaporation 
processes within the ISM. 

As an example we show the dynamical and chemical evolution of a star forming 
dwarf galaxy with a total baryonic mass of 2 • 10 ®Mq. After a moderate collapse 
phase the stars and the molecular clouds follow an exponential radial distribution, 
whereas the diffuse gas shows a central depression as a result of stellar feedback. 
The metallicities of the galactic components behave quite differently with respect 
to their temporal evolution as well as their radial distribution. Especially, the ISM 
is at no stage well mixed. 
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1. The 3d chemodynamical code and first results 


Since several years smoothed particle hydrodynamics (SPH) calcula¬ 
tions have been applied successfully to study the formation and evo¬ 
lution of galaxies. Its Lagrangian nature as well as its easy implemen¬ 
tation together with standard N-body codes allows for a simultaneous 
description of complex dark matter-gas-stellar systems (e.g. Navarro &: 
White, 1993; Mihos & Hernquist, 1996). Here we present simulations 
based on our 3d chemodynamical code. This code includes many com¬ 
plex effects such as a multi-phase ISM, cloud-cloud collisions, a drag 
force between different ISM components, condensation and evaporation 
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Figure 1. The radial distribution of the cumulative mass (left) and the surface 
density (right) for the different components in the central region of the model galaxy 
after 1 Gyr. 


of clouds (CE), star formation (SF) and a stellar feedback (FB). This 
code is a further development of our single phase galactic evolutionary 
program (Berczik, 1999; Berczik, 2000) including now different gaseous 
phases. 

In our new (multi-phase gas) code we use a two component gas 
description of the ISM (cold “clouds” and “smooth” hot SPH phase) 
(Theis et ah, 1992; Samland et ah, 1997). The basic idea is to add a cold 
cloudy component to the smooth and hot gas (10"^ - 10^ K) described 
by SPH. The cold clumps are modeled as N-body particles with some 
“viscosity” (Theis & Hensler, 1993) (cloud-cloud collisions and drag 
force between clouds and hot gas component). The cloudy component 
interacts with the surrounding hot gas also via condensation and evap¬ 
oration processes (Cowie et ah, 1981; Koppen et ah, 1998). In the code 
we introduce also star formation. The “stellar” particles are treated 
as a dynamically separate N-body component. The cloud component 
forms stars. These stars return chemically enriched gas material and 
energy to both gaseous phases. 

As a test of our new code, we calculate the evolution of an isolated 
star forming dwarf galaxy. The initial total gas content of our dwarf 
galaxy is 2-10® Mq (80 % “COLD” -h 20 % “HOT”) which is placed in¬ 
side a fixed dark matter halo with parameters rg = 2 kpc and po = 0.075 
Mq/pc^, using density profile of Burkert, 1995. With these parameters 
the dark matter mass inside the initial distribution of gas (20 kpc) 
is ~ 2 • 10^^ Mq. For the initial gas distribution we use a Plummer- 
Kuzmin disk with parameters a = 0.1 kpc and b = 2 kpc (Miyamoto & 
Nagai, 1975). The gas initially rotates in centrifugal equilibrium around 
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Figure 2. The temporal evolution of the mass (left) and mass exchange rate (right) 
for the different components of the model galaxy. 


the z-axis. We choose the dwarf galaxy as an appropriate object for 
our code, because in this case even with a relatively “small” number 
of cold “clouds” (~ 10 ^) we achieve the required physical resolution 
for a realistic description of individual molecular clouds ( 10 ® Mq) 

as a separate “COLD” particle. In the simulation we use Nhot = 10'^ 
SPH and Ncom = 10^ “COLD” particles. After 1 Gyr more then 10'^ 
additional stellar particles are created. 

In Fig. 1 we present the mass and surface density distribution of the 
different components in the central region of the model after 1 Gyr of 
evolution. In the region up to ~ 2 kpc the baryonic matter dominates 
over the DM. The surface density of the stars can be well approximated 
by an exponential disk with a scale length of 0.55 kpc. In the distribu¬ 
tion of hot gas, we see a central “hole” 1 kpc), as a result of gas 
blow-out from the center mainly due to SNII explosions. 

In Fig. 2 we present the evolution of the mass and the mass exchange 
rate of the different components. The SFR (i.e. dMsTAR/dt) peaks to 
a value of 1 M 0 yr“^ after 200 Myrs. Afterwards it drops down to 0.2 
M 0 yr“^ within several hundred Myrs. Another interesting feature is 
the behaviour of the hot gas phase mass exchange. After the initial 
violent phase of condensation an equilibrium is established which gives 
a hot gas fraction of about 10 % of the total gas mass. 

The metal content of the diffuse gas and the clouds differs signifi¬ 
cantly over the whole integration time (Fig. 3). Due to SNII and SNIa 
events the metallicity of the hot phase exceeds that of the clouds by 
almost one order of magnitude. The clouds mainly get their metals by 
condensation of the hot phase. The central metallicity plateau (up to 
1 kpc) of the cold component is explained by the fact, that condensation 
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Figure 3. Temporal evolution of the metallicities (left) and their radial distribution 
after 1 Gyr (right). Individual metallicities of newly born stars are marked by dots. 


is not very efficient in that region, e.g. because of the central “hole” of 
the very hot diffuse ISM and a lack of metal enriched material there. 
Moreover, the conditions in the center lead mainly to evaporation of 
clouds which also prevents the mixing with the metal enriched hot gas. 
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